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ABSTRACT 


Methods are described for calculating significant factors 
in the observational environment of orbiting astronomical 
satellites. These factors must be considered in the pro- 
cess of scheduling observations and in data reduction. 
Subroutines which perform these calculations are 
described. 
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THE OBSERVATIONAL ENVIRONMENT OF ASTRONOMICAL SATELLITES 
AND RELATED SOFTWARE SUBROUTINES 

Section 1 
INTRODUCTION 

Astronomical observations made from oribiting satellites are the end pro- 
duct of high cost space astronomy programs. To be cost effective these pro- 
grams require that the scientific operation of the observatory must be highly 
efficient and that the resulting data be of highest possible quality, free from the 
contaminating background effects which are ever present in the orbital observa- 
tional environment. 

The data contaminating backgrounds which most seriously constrain OAO 
observations are: 

1. Background radiation due to the concentration of high energy protons 
and electrons in the earth bound region known as the South Atlantic Anomaly. 

The OAO orbit passes directly through this region several times per day. Both 
OAO-II and OAO-C are seriously contrained by this background loosing about 
30% of their total available observing time. 

2. Background light which is scattered off the OAO light baffles and 
enters the telescope optics whenever the spacecraft pointing axis comes close 
to bright sources such as the sun the moon and the sunlit earth. These 
scattered light backgrounds consume an additional 30% of the available ob- 
serving time on OAO. 

These background constraints leave only about 30% of the time when the 
orbital environment is reasonably "clean" for the acquisition of reliable data. 

The orbital observational environment is of course a highly dependent 
function of time. It is a compound function of the relative positions of the sun, 
the moon, the earth, and orbit of satellite and the position of the satellite in 
its orbit. 

In addition there are other constraints called operational restrictions. 

These are due to spacecraft thermal and power limitations which prohibit pointing 
the spacecraft within certain critical regions. These constraints are a function 
of the position of the sun and the precession of the orbit. 

In order to successfully carry out the scientific objectives it is essential 
that the scheduling of astronomical observations be done efficiently and with care 
to insure that maximum use is made of the available viewing time and that the 
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data acquired during this time is of highest possible purity. The scheduling of 
observations and the subsequent reduction of data also must maintain a pace 
equal to or greater than the maximum observational rates which are controlled 
by the orbital environmental constraints and the spacecraft operational 
restrictions. 

It is clear that computer software systems must be used to accomplish 
the above stated requirements. Orbital elements, spacecraft parameters and 
observational requirements are input to the system. The subroutines described 
in this report compute the orbital environment as a function of time. 

By a sequence of calls to these subroutines the schedule driver assemblies 
the critical beginning and end times for periods of environmental contamination 
and constructs a timeline. Other programs then integrate this environmental 
timeline to accomplish the various tasks needed to schedule observations during 
periods which conform to the data requirements of a particular observation. 

In addition to its value in scheduling orbital observations the environmental 
timeline is essential to the comprehensive reduction of data. The timeline is 
an accurate concise minute by minute history of conditions under which the data 
was acquired. It provides the necessary information required to evaluate the 
quality or the contamination of the data and it rapidly facilitates the subsequent 
background corrections. The most effective way to correlate the environment 
with the data is to input the timeline to the data reduction program. Examples 
of environmental timeline formats used in OAO scheduling and data reduction 
programs are shown in Figures 1 and 2. 


Section 2 

METHODS OF SOLUTION 


2. 0 General Information 

Described in this section are the methods employed in the calculation of 
the significant factors of the OAO environment. The expressions which follow 
define the geometry of the orbit both in celestial and geocentric frames, and 
consequently permit expansion to include factors which have not been considered 
significant in the past but which, with increasingly sophisticated instruments, 
may become significant in the near future. 

While these methods were developed specifically for the OAO-B space- 
craft it has been demonstrated that they are applicable to any circular orbit 


2 



and have been used to generate a simulated orbit for the proposed IUE satellite 
which will have a synchronous altitude. 

2. 1 Position of the Principle Bodies 

The positions of the Sun and Moon with respect to the satellite are assumed 
to be the same as with respect to the Earth. Ephemeris equations that generate 
Right Ascension and Declination for the Sun and Moon were developed for use in 
the OAO-B Scheduling System (Reference 1). These equations are functions of 
time expressed in Julian Day minus 2.438 million. Unit vectors are calculated 
from R. A. and Dec. using the relationships: 


where 


U J 

= cos 

a. 

cos 

6 

u 2 

= sin 

a 

cos 

6 

U 3 

= sin 

6 




a is Right ascension 
6 is Declination 


2. 2 Calculation of Qccultation Period 

The criteria for determining target occultation by the Earth are the 
same as those used in the ground support programs to command the satellite. 
Rather than using the true Earth size, a safety region is created by increasing 
the true angular radius of the Earth by an amount specified as an input. Usually 
this is about 14 degrees. The total angular radius of the Earth is then: 



where 


Re is the Earth's mean radius 
A is the semi-major axis of the OAO orbit 
D is the width of the safety region in radians 

The calculation of occultation assumes a spherical Earth and a circular 

orbit. 


* 
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First, we calculate the time, t eo , corresponding to the nodal passage of 
the orbit for which the elements are given. 


teo — to — 


We 


+ M,, 


where 

t e is the epoch in days 
o) e is the argument of perigee at epoch 
M e is the mean anomaly at epoch 
p angular motion in radians per day 

Then, the number of orbits completed since epoch, an integer, may be found by, 


t is the present time 

„ 2 7T .... 

P = -p- period in days 

The right ascension of the ascending node is then 
0 - « e + Nd 

where 

S2 e is the right ascension of the ascending node at epoch 

d is the precession of the node in radians per orbit 

Next we must find the azimuth of the vector along the negative boresight 
axis (the vector from the target to the satellite) projected into the orbital plane. 

The x and y direction cosines of the boresight vector in the orbital plane 
are given by 


a = cos 6 1 cos (n-at ) 

and 

b = cos i cos 6 t sin ( J2 - Q! t ) 
+ sin i sin 6 t 
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where 


and 


by 


oi t is the right ascension of the target 
S t is the declination of the target 

i is the inclination of the orbit. 

The azimuth of the negative bores ight vector is then: 



The points of emersion and immersion differ from e by an angle rj given 


r] = cos' 1 


‘ cos y ' 
Va 2 + b 2 . 


where 


y, a and b are defined above. 

Thus we have the anomaly at emersion of the target: 

0 r = € + TJ 


and at immersion 

0, = € - f? 

The times corresponding to and /3 S are obtained by: 
t r = teo + NP + i3 r /p 
ts = teo + NP + 0 s /p 

where all terms have been previously defined. 

A complete mathematical treatment of this problem may be found in 
Reference 2 from which this solution was derived. 
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2. 3 Calculation of Radiation Level in the South Atlantic Anomaly 


When the satellite passes through the South Atlantic Anomaly, the com- 
paratively high-level of background radiation may have a significant effect on 
the data. For this reason, it would assist the astronomer in analyzing the data, 
to know approximately what the background flux intensity was when the data was 
taken. Reference 3 contains world maps of the South Atlantic anomaly derived 
emprically and displayed parametrically with respect to altitude and threshold 
energy. For this application, the map of proton flux above 15 MEV for an 
altitude of 800 KM was used. This map gives three contours corresponding to 
proton fluxes of 10, 100 and 1000 protons per cm 2 per sec. Each of the three 
contours was fit with a polygon. If the latitude and longitude of the sub-satellite 
point are known, it may be determined whether the point lies inside any of the 
contours. 

Determination of the sub-satellite point begins with the calculation of 
satellite position in celestial coordinates. The right ascension of the ascending 
node is given by: 


= fi e + (t — t e ) d 

where 

S2 e is the right ascension of the ascending node at epoch. 
t e is the epoch of the orbital elements in days 
t is present time in days 

and d is precession of the node in radians/day. 

The azimuth of the satellite position measured from the ascending node is: 
A =co e +M e + (t-te)p 

where 

w e is the argument of perigee at epoch 
M e is the mean anomaly at epoch 
pis the angular motion in radians per day 
t and t e as above. 
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To calculate the longitude, we make use of an ephemerls equation which yields 
Greenwich Siderial Time as a function of Greenwich Mean Time, or: 

t s = f (t) 

Since GST is by definition the right ascension of the Greenwich Meridian, the 
satellite longitude, 9 , is: 

9 = a - t s 

The sub-satellite point (9, 0) may now be compared with the contours to find 
the approximate proton flux. The solution was derived from work done on the 
OAO-B Scheduling System (Reference 1). 

2. 4 Calculation of Percent of Visible Earth Illuminated by Sunlight 

It is useful to the astronomer analyzing the data to know what percent of 
the Earth which is visible from the spacecraft is illuminated by the sun. Sub- 
routine PELS which performs this calculation is available from the OAO Space- 
craft Control Programming Section (NASA Code 542) and the mathematical 
analysis is described in detail in Reference 4. This routine was developed for 
use in the OAO-B Scheduling System (Reference 1). 
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2. 5 Calculation of Angular Distance to Nearest Illuminated Point of Earth 

In order to evaluate the effects of an illuminated Earth on the data, the 
astronomer must know the angle between the experiment bo resight and the 
illumination source. Since the Earth subtends over one-fourth of the total 
solid angle as seen from the spacecraft, the most meaningful quantity to calcu- 
late would be the nearest lit point. This calculation is performed by Subroutine 
DTSE which is available from OAO Spacecraft Control Programming Section 
(NASA Code 542). The documentation for this subroutine was not yet available 
at the time of writing, however, the analysis may be obtained from the same 
source as Reference 4. 


Section 3 

DESCRIPTION OF SUBROUTINES 
3. 0 General Information 


The four routines described in this section form a package which ma y 
be incorporated into the driver of scheduling or data reduction system to ac- 
curately simulate the OAO environment as a function of time. 

It is important to preserve the order of execution of these routines be- 
cause certain intermediate quantities are generated in one routine and passed 
through the program interface to another. In general OCCULT should be 
executed first to read in orbital data and then for each time element of the 
simulation OCCULT should be executed first, then ANMLY followed by PELS 
and DTSE the latter two which may interchanged. 

3. 1 Subroutine ANMLY Functional Description 

Subroutine ANMLY generates an estimate of the radiation level experienced 
by the satellite as it passes through the South Atlantic anomaly. This is done 
by calculating the latitude and longitude of the sub-satellite point using the 
method described in Paragraph 2. 3 and then testing to see whether the point 
lies within any of three map contours of the anomaly area. 

The map contours represent boundaries of regions where proton radiation 
exceeds 10, 100 and 1000 protons/cm 2 /sec. These particular levels were se- 
lected for OAO-B and can be easily re-defined for any other satellite. 
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A single integer is set to zero if the sub-satellite point is not within any 
of the contours, and is set to 1, 2 or 3 if the point lies within the 10, 100 or 
1000 proton/cm 2 /sec. contour respectively. This integer then is approximately 
equal to log 1Q F where F is the proton flux. 

A call to subroutine OCCULT must be executed prior to calling ANMLY 
because OCCULT reads the orbital elements from a disk dataset and calculates 
several quantities used by ANMLY. 

Also, ANMLY must be executed prior to calling PELS or DTSE for each 
simulation time element because ANMLY generates the satellite position 
direction cosines used by those subroutines. 

Subroutine ANMLY may be easily modified to estimate cosmic ray flux 
as a function of magnetic latitude which can be generated from the latitude and 
longitude of the sub-satellite point. 

Program Interface 


Calling sequence: 
CALL ANMLY (AJID, NOM) 


AJID 

TO 

EJD 

DEGDAY 

ANODE 

ANGE 

TPP 

FINC 


Input Items: 

Julian Day of digital group being processed. 

(common) ratio of earth radius to orbital semi-major axis. 

(common) Julian Day of epoch of orbital elements. 

(common) orbital precession of node in degrees/day. 

(common) Right Ascension of ascending node at epoch. 

(common) sum of the mean anomaly at epoch and the argument 
of perigee at epoch. 

(common) angular motion in radians/day. 

(common) orbital inclination in radians. 
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Output Items: 


NOM approximate integer value of log base 10 of proton flux. 

PSAT (common) direction cosines of satellite position in inertial 

coordinates. 

Common Blocks Referenced: 


COMOCL 

3. 2 Subroutine OCCULT Functional Description 


Subroutine OCCULT is a modified version of Subroutine OCTIM of the 
OAO-B Scheduling System (Reference 1). Subroutine OCCULT calculates the 
positions of the sun and moon and computes the times of emersion and immersion 
of occultatidhs of the target by the earth (see Paragraphs 2. 1 and 2. 2). 

It is interesting to note that the capability readily exists for inclusion of 
a correction for wavelength aberration due to motion of the spacecraft normal 
to pointing axis since all of the orbital parameters are present. 

The first call to OCCULT reads the orbital elements dataset and com- 
putes the earth's angular radius, the Julian Day of Epoch and the orbital period. 

Subsequent calls to OCCULT return the rise and set times of the target 
for a specified orbit number and target position. An alternate entry point, 
OCCTIM, is provided which is entered with Julian Day rather than orbit num- 
ber, the latter being calculated upon entry. 

Entry point OCNOD returns the time of nodal passage in Julian Days for 
a specified orbit number. 

The OCCULT and OCCTIM entries have an intermediate results print 
flag in the argument list which if greater than zero, causes a minute-by-minute 
timeline plot of the orbit showing the unocculted intervals by the letter 'U' and 
the occulted portions with a dash. 
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Program Interface 


Calling sequence: 

CALL OCCULT (NORB, RISE, SET, RAS, DECS, LMN, ANAME) 
CALL OCCTIM (AJID, RISE, SET, RAS, DECS, LMN, ANAME) 
CALL OC NOD (NORB, TIM) 



Input Items: 

NORB 

orbit number. 

RAS 

Right Ascension of target in radians. 

DECS 

declination of target in radians. 

LMN 

intermediate printout flag. 

ANAME 

starname 

AJID 

current time in Julian Days. 


Orbital elements and principal body sizes read from permanent disk 
dataset (see Section on Program Inputs). 

Output items: 

Julian Day of target emersion from occultation by earth. 
Julian Day of target immersion behind earth. 

(common) Julian Day of nodnal passage. 

(common) angular motion in radians per day. 

(common) Julian Day of Epoch. 

USUN (common) direction cosines of sun position. 


RISE 

SET 

ANJD 

TPP 

EJD 
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Subroutine called: 


JDAY 

TIME (see DATE) 

Common blocks referenced: 

COMOCL 

3 . 3 Subroutine PELS Functional Descript ion 


Subroutine PELS was developed for OAO-B with the help of the Space- 
craft Control Programming Section NASA Code 514 (see Reference 4). PELS 
calculates the percentage of the earth, as seen from the spacecraft, which is 
lit by the sun. 


Program Interface 


Calling sequence: 

CALL PELS (USUN, PSAT, A, B, TOL, IPR) 

Input Items: 

USUN unit vector to sun in inertial coordinates. 


PSAT position vector of satellite with respect to the earth in units of 

earth radii. 


graying tolerance in radians used to handle quadrature cases 
where circular functions are badly behaved, for this applica- 
tion set to . 001. 


intermediate printout flag; 1-no printout, 0-print out arguments 
only, -1- full debug printout. 

Output items: 

A solid angle in steradians subtended at satellite by sunlit earth. 

B A expressed as percentage of solid angle of full earth. 

Subroutines called: 

DQSF (Scientific Subroutine Package) 
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3.4 Subroutine DTSE Functional Description 

Subroutine DTSE was developed for OAO-B with the help of the Spacecraft 
Control Programming Section NASA Code 514. DTSE calculates the angular 
distance from the spacecraft boresight axis to the nearest sunlit point of earth. 

Program Interface 

Calling sequence: 

CALL DTSE (USUN, PSAT, ATUDE, DST, IPR) 

Input items: 

USUN unit vector to sun in inertial coordinates. 

PSAT position vector of satellite with respect to earth in units of 

earth radii. 

ATUDE unit vector of spacecraft boresight in inertial coordinates. 

IPR intermediate printout flag; -1 no printout; 0 arguments only; 

1 partial debug printout; 2 full debug printout. 

Output items: 

DST angular distance in radians from boresight to nearest sunlit 

point of earth. 

Subroutines called: 

CROSS SINSLV 

DOT F 

NORM DRTMI (SSP) 

3. 5 FUNCTION JDAY Functional Description 

Function JDAY returns integer Julian Day from the month, day of month 
and year or from the day of year and year. 
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Program Interface 


Calling Sequences: 

J = JDAY (ID, IM, IYR) or J = JDAY (IDY, IONE, IYR) 

Input Items: 

ID, IM, IYR day of month, month, and year minus 1900. 

IDY day of year. 

IONE integer variable containing the integer value 1. 

Output items: 

J integer variable into which integer Julian Day is to be stored. 

3. 6 Subroutine TIME Functional Description 

Subroutine TIME converts floating point Julian Day into calendar date and 
time in hours, minutes and seconds. Day of year is also calculated for con- 
venience. 

Program Interface 
Calling Sequence: 

CALL TIME (AJID, ID, IM, IY, IGDAY, IHR, IMIN, ISEC) 

Input Items: 

AJID Julian Day to be converted 

Output Items: 

ID, IM, IY, IGDAY, IHR, IMIN, ISEC - day of month, month, year minus 1900, 
day of year, hour, minute and second corresponding to Julian Day AJID. 
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Figure 1. Example of Utilization of Environmental Data in 
Observation Scheduling 
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Figure 2. Example of Utilization of Environmental Data in Data Reduction 



